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Abstract 

The hysteresis loop in the zero-temperature random-field Ising model ex- 
hibits a critical point as the width of the disorder increases. Above six dimen- 
sions, the critical exponents of this transition, where the "infinite avalanche" 
first disappears, are described by mean-field theory. We expand the critical 
exponents about mean-field theory, in 6 — e dimensions, to first order in e. 
Despite e = 3, the values obtained agree reasonably well with the numerical 

values in three dimensions. 
PACS numbers: 75.60.Ej, 64.60.Ak, 81.30.Kf 



1 



Typeset using REVTgX 



In a previous paper [0, we modeled hysteresis in magnetic and martensitic systems 
using the random-field Ising model at zero temperature. The model exhibited two features 
characteristic of these systems: the return-point memory effect and avalanche-generated 
noise. (The noise is called Barkhausen noise in magnetic systems and acoustic emission in 
martensites.) We also discovered a critical point, separating smooth hysteresis loops at large 
disorder where all avalanches are finite, from discontinuous hysteresis loops at small disorder 
where one avalanche turns over a fraction of the whole system. 

Here we study this critical point in an expansion about mean-field theory. Figure [l|(a) 
shows a schematic of the phase diagram for our model defined by eq. @ below. The vertical 
axis H is the external field. The horizontal axis R is the width of the probability distribution 
of the random fields fi acting on each spin. The bold line represents the location Hc{R) at 
which the infinite avalanche occurs, when the field H{t) is adiabatically increasing from an 
initial state where all spins were pointing down. At small disorder, the first spin to fiip easily 
pushes over its neighbors, and the transition happens in one burst (the infinite avalanche). 
At large enough disorder, the coupling between spins becomes negligible, and most spins 
fiip by themselves: no infinite avalanche occurs. At a special value of the randomness 
R = Rc the infinite avalanche disappears. We find a critical point with two relevant variables 
r = {Rc — R)/Rc and h = {H — Hc{Rc)) lU- At this point we find a universal scaling law 
for the magnetization m = (M — Mc{Rc)) 

m ~ |r|/^A^±(VkP^), (1) 

where the ± refers to the sign of r. 

We use a soft-spin version of the random field Ising model, whose energy at a given spin 
configuration {sj} is 

= - ^ JijSiSj - fiSi + Hsi - V{si) , (2) 

ij i 

with the linear cusp-potential V{si) defined 0] through 
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k/2 {si + 1)2 for s < 
k/2 {si - 1)2 for s > 

Here, k > is the local curvature of the potential. The spins are coupled ferromagnetically 
by a nearest neighbor interaction Jij = J/ z > 0, z being the coordination number of the 
lattice. We demand that A;/J > 1 to ensure stability of the system. H is a homogeneous 
external magnetic field; the are randomly chosen from a Gaussian distribution p(/) of 
standard deviation R. We study this system at zero temperature. It turns out for given H 
and fi that there are many metastable states; which one of these the system picks depends 
entirely on its history (i.e. the way the external magnetic field H was varied at earlier times). 
We will study the history of a monotonically but adiabatically increasing external magnetic 
field. We impose purely relaxational dynamics, as defined by the equation of motion 

where we have absorbed the friction constant into the definition of the time t. 

a. Formalism: We use the formalism of Martin, Siggia and Rose to rewrite the 
problem as a path integral for a generating functional Z, and then expand this functional 
about mean field theory. This is done in analogy with the calculation for CDWs by Narayan 
and Fisher We impose the dynamical equation on the path integral at each time 

t by introducing it as a 5-function constraint using the well-known identity 27r5(/(s)) = 

Z = J[ds][ds]J[s]exp{S) (4) 



where 



S = ^J dtJ2s,itMs,it)-J2J^eSe-H-fj + ^) (5) 

j i J 

Derivatives of Z can be used to calculate the response functions and dynamic correlation 
functions for our model. J[s] is a functional Jacobian, chosen such that Z is unity, indepen- 
dent of the fj. This allows us to average over the disorder without fancy tricks (like replica 
theory) . 



We choose a particular regularization for the time integral. The simplest choice |^ is 
to require a force at time t to have an effect only after some time 6t. That leaves us with 
J[s] = 1. We now do an average over the random fields /«, denoted by ( )/, leading to the 
averaged generating functional 

Z = J[ds][ds]{expiS))f. (6) 

To expand about mean-field theory, we need change variables from sj and sj to the local 
fields rij = (1/ J) J2e JjeSe at the sites (fluctuations about whose mean values we shall study). 
We do so by introducing another auxiliary field "fjj, and absorb a factor i in its definition, so 

where Zj[rij,'fjj] is a local functional 

h ^^j]= J [^^] [^^] (e^P •S'j ) / , (8) 

and 

We will now expand about the mean-field solution tjq (the local field configuration about 
which the log of the integrand in equation (|^ is stationary). Shifting the definition of 1] to 
rj — rjo so that {tj) j = leaves one with the generating functional 



J[d7^mexp{S) (10) 



with an effective action 



r °° 1 r 

j,l J m,n=0 

'^mn('^l) ■■■:tm+n)Vj(tl) ' ' ' Vji.^m)Vji.^m+l) ' ' ' Vji.^m+n) ■ 

Here, the Umn are the derivatives of In Zi with respect to the fields fjj and rij, and thus are 
equal to the local, connected responses and correlations in mean-field theory: 



Local (/) means that we do not vary the local field (?7o)j in the mean-field equation 

dtsj{t) = JivoUt) + H + f,- + Je{t) (13) 

when we perturb with the infinitesimal force Je{t). The force Je{t) is only allowed to increase 
with time, consistent with the history we have chosen. (For example for mi i(t,t') we add a 
force Je 9(t — t,) in eq. (p!3|), with 0(t — t/) being the step- function, and take the derivative 
of {s(t))f^i^c with respect to e and t,.) 

b. RG treatment: We consider the 1777 term in the action (involving J^^^ and Ui 1) as the 
propagator in the RG treatment. Now we take some long- wavelength and low-frequency 
limits in analogy to [|],||. For small wave vectors J~^{q) ~ (1/^) + <^2Q'^5 and we rescale to 
give JJ2 = 1. We take the low-frequency part of the propagator, by Fourier transforming 
the fjri term in time, expanding to first order in uo and Fourier transforming back. The 
propagator (the 7)7] term in the action) is thus (after rescaling) 

- jd^j dtfi{-q,t)[-dt + q'-x-'/JHq,t). (14) 

The bare value of x is the static response, calculated in mean-field theory, to a monotonically 
increasing external magnetic field 

X = 1/ (2J^p{-JM - H + k) ~ J{k - J)/k) , (15) 

where M is the magnetization at the external magnetic field H. 

We use the Wilson-Fisher renormalization group transformation: In each step we inte- 
grate out modes of all frequencies and wave vectors within an infinitesimal wave vector shell 
0. We rescale through x = bx', t = b^t'. We choose the rescaling of the fields such that the 
term of the propagator and the ^2,0 term remain unchanged, since to first order in e they 
have no loop corrections. Thus fj = b'i'^ff' and r] = b^i^'^rj' [^]. Without loop corrections 
this implies z = 2. Keeping in mind that the d/dt{t) in eq. ([12]) rescale like we arrive at 



^mn ~ To lowest Older in e = dc — d, the only relevant terms are those 

which do not flow to zero under rescaling at the upper critical dimension d^ (we will see that 
dc = Q for our critical point), ui^ is trivially zero because we expand around the stationary- 
point. The Ui^2 term, in the static limit, has bare value w = —2J^p'{—JM — H + k): it 
becomes relevant for d < 8. The Mi_3 term starts at m = 2J^p"{—JM — H + k) in the static 
limit, and is relevant for d < 6. Finally, the M2,o term stays marginal. In the static limit 
that we consider M2,o can be treated as a constant. We are left with the effective action 

^ = -|A/ dtv{-q,t)[-dt + q'-X-'/J]viq,t) + J2[i^/'^) J dtf^,{t){r^^{t)fw 

i 

+ (1/6) / dtfjj{t){'n,{t)fu + {1/2) [ dti f dt2flj{ti)fl,{t2)u2,o{tut2)] (16) 



Our e-expansion can be applied not only to the critical point {Rc, Hc{Rc)), but to the 
entire line Hc{R) at which the inflnite avalanche occurs. In mean-fleld theory, the approach 
to this line is continuous, with a power-law divergence of the susceptibility x ^i-nd precursor 
avalanches of all scales. Above 8 dimensions the action is purely quadratic at the flxed 
point, and the inflnite avalanche line (where 1/x = and w = — 2J^p'(— JM — H + k) 7^ 0) 
presumably remains critical. For d = 8 — e, flgure |^(a) shows the correction to vertex w to 
flrst order in e. The incoming lines at a vertex stand for rj operators, and the outgoing lines 
are fj operators. The the low-frequency form of the propagator is approximately 6{t — t')) ^ 
but we have to observe causality; an example of a diagram forbidden by causality is given 
in flgure ^(b). Applying the usual approximations 0, we obtain for the recursion relation 
for w to 0(e): 

w'/2 = i^w/2 + {u2,o/2){w/2Y8/{47r)' J\ql/{q' - x"V^)'} (17) 

Writing b^~'^/'^+^^ = ll^'^l'^) = 1 + e/2 log 6 and performing the integral over the momentum 
shell A/b < q < A leaves us with the recursion relation: 

w'/2 = w/2 + {w/2) (e/2 + M2,o(w/2)'4/(47r)^ log h) (18) 

Since ^2,0 > this means that for e > there are only two flxed points with w' = w: either 
w = 0, which we will discuss in the next paragraph, or w = 00. We see that under the 
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recursion relation (|T8]) any system that has a bare value w ^ when 1/x = will flow 
to the fixed point w = (yo. We interpret this as indication that the transition is a first 
order transition for d < 8. Indeed, in three dimensions the simulation showed a first-order 
transition without critical fiuctuations for these systems. 

The critical point we are interested in here is the fixed point where w = 0. At d = 6 the 
first non-quadratic contribution u becomes relevant, i.e. the upper critical dimension [0] for 
the critical endpoint is 6. We now compute, to 0(e), the corrections to the recursion rela- 
tions. The relevant diagrams are shown in figure |^(c). Figure |l|(b) shows the corresponding 
RG fiows in the plane. 

The loop-corrections look very similar to the loop-corrections in the usual Ising model 
in d — 2 dimensions. In fact, to 0(e), they are the same. This can be seen either by direct 
computation (see caption figure (0)), or by noticing that the Feynman rules for our diagrams 
are the same as those for the finite-temperature random-field Ising model (except that 
we have extra vertices which are irrelevant to 0(e)). This latter model has been mapped to 
all orders in e onto the regular Ising model, using supersymmetry and other arguments. This 
analogy tells us that to 0(e) we get the same RG fiows (figure 1) and the same corrections 
to our exponents as one finds in the usual Ising model in — 2 dimensions (see table). 

This mapping does not extend to the next (e^) term in the series: figure H(d) shows 
a correction of O(e^) to the vertex ^2,2, which then contributes in O(e^) to the propoga- 
tor. This is comforting, as otherwise the critical properties of our model in c? = 3 would 
have completely mapped onto the d = 1 thermal (non-random) Ising model, which has no 
finite-temperature phase transition at all. Indeed, this was a substantive concern for the 
thermal random-field Ising model, which despite the correspondence above was proven to 
have a transition in = 3: the e-expansion for that model summed over physically incorrect 
metastable states. By controlling the history of the external field (as in PJ§]), we have been 
careful to specify the particular metastable state in our calculations. 

The e-expansion for our model is technically much simpler than that for other disordered 
extended dynamical systems: e.g. interface or charge- density wave P| depinning, where 
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an infinite family of relevant operators made necessary a functional renormalization group. 
The relative ease of our calculation may make possible further extensions: calculating the 
corrections to the equations of state, calculating the history-dependent critical behavior, or 
addressing the avalanche distributions. 

We would hke to thank O. Narayan and D. S. Fisher for advice and consultation, and 
S. Kartha, J. A. Krumhansl, M. E. J. Newman, B. W. Roberts, S. Ramakrishna, J. D. Shore, 
and J. von Delft for helpful conversations, and NORDITA where this project was started. 
We acknowledge the support of DOE Grant #DE-FG02-88-ER45364. 
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TABLES 



exponents 


e expansion 
with e = (i — 6, ate = 3 


simulation Q 
in 3 dimensions 




2 - e/3 = 1 


1.0 ±0.1 




0.5 - e/6 = 


0.17 ±0.07 


[55 


1.5 + 0(e2) = 1.5 


2.0 ±0.3 


6 


3 + e = 6 


(around 12) 


TABLE I. Universal exponents for critical behavior in hysteresis loops. The exponents 
(3 and 6 tell how the magnetization scales with r and h, respectively, equation (^). v is the 
correlation length exponent, measured (numerically) using finite-size scaling. 
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FIGURES 

FIG. 1. Phase diagram and flows (schematic), (la) The vertical axis is the external 
field H, responsible for pulling the system from down to up. The horizontal axis is the width of 
the random-field distribution R. The bold line is Hc{R), the location of the infinite avalanche 
(assuming an initial condition with all spins down and a slowly increasing external field). The 
critical point we study is the end point of the infinite avalanche line {Rc, Hc{Rc))- 

Using the analogy with the Ising model (see text) we also show the RG fiows around the 
critical point. Here we ignore the RG motion of the critical point itself: equivalently, the figure 
can represent a section through the critical fixed point tangent to the two unstable eigenvectors 
(labeled h and r). Two systems on the same RG trajectory (dashed thin lines) have the same 
long- wavelength properties (correlation functions ...) except for an overall change in length scale, 
leading to the scaling collapse of equation (jl]). The r eigendirection to the left extends along the 
infinite avalanche line; to the right, we speculate that it lies along the percolation threshold for up 
spins (see reference |^). 

(lb) 0(e) RG fiows below 6 dimensions in the (x^^,'u) plane (see text). Linearization around 
the Wilson-Fisher (WF) fixed point yields the exponents given to 0(e) in the table. In the vicinity 
of the repulsive u = = (MFT) fixed points one obtains the old mean-field exponents. 

FIG. 2. Feynman diagrams. The perturbative expansion about mean-field theory is pre- 
sented here by Feynman diagrams, (a) The correction to 0(e) to the vertex w in an expansion 
about 8 dimensions, see eq. ( p!^ ) in the text, (b) An example of a diagram forbidden by causal- 
ity, (c) The relevant corrections to first order in e = 6 — d for the constant part x^^ / J the 
propagator and for u. Using the same techniques that lead to eq. (^) we find the following re- 
cursion relations: (x^'^/J)' = &^(x"V-'^ + ^2,oV(47r)^A^(l - l/b^)/4: + U2,ou/{4tt)''^{x''^ / J)logbJ 
and u' = u + u(^e + 3/(47r)^ti2,o'") log 6. U2fl does not get any loop corrections of 0(e). (d) An 
example of a correction to the vertex ^2,2 which contributes only to O(e^), which is not present in 
the regular Ising model (or the thermal random field Ising model). 
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